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ABSTRACT 



In this thesis, the algorithms for an A-7E aircraft perfor- 
mance calculator were developed and then implemented on three 
small data processors of different programming levels and storage 
capabilities . 

The utility of data is a function of several variables in- 
cluding accuracy and availability. The problem of retrieving 
performance data from the Naval Air Training and Operating Pro- 
cedures Standardization (NATOPS) Manuals is significantly les- 
sened by the devices demonstrated in this investigation. Nine 
performance chart groups, yielding data usually considered 
necessary for flight, were reduced to a series of analytical 
expressions. These analytical expressions were demonstrated to 
reproduce NATOPS Manual data to a high degree of accuracy. 

Implementation was demonstrated on a desk computer, a hand 
held calculator and a microprocessor. 
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I . INTRODUCTION 



The Naval Air Training and Operating Standardization 
(NATOPS) Manual is the official standard of the United States 
Navy for "...information on all aircraft systems, performance 
data, and operating procedures required for safe and effective 
operations." [1] 

The purpose of this thesis was to develop algorithms of 
the more often used NATOPS performance charts for the A-7E 
aircraft, examine their accuracy and implement them on small 
data processors that might be adaptable to shipboard or aircraft 
onboard use. The interpretation of NATOPS performance charts 
is an error prone and time consuming procedure even for experi- 
enced users. The need for a system to eliminate this laborious 
process has been fully documented in a thesis completed in June 
1978 by LCDR W.M. Siegel [2]. In his investigation, LCDR Siegel 
devised an efficient procedure to develop algorithms from the 
NATOPS performance charts and exercised this procedure on the 
problems of "Takeoff Ground Roll Distance" and "Takeoff Airspeed". 

This investigation is an extension of the aforementioned 
work. The original scope of this investigation was to develop 
algorithms for eleven of the most often used performance problem 
chart groups and implement them on the Texas Instruments- 5 9 
(TI-59) hand held calculator (.HHC) . All of the NATOPS perfor- 
mance charts were not reduced because of research time 
limitations. Of the eleven performance chart groups studied. 
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two performance problems, "Time to Climb" and "Fuel Required 
to Climb" were rejected because of implementation difficulties 
on the TI-59 HHC Cdiscussed fully in "Development Difficulties") . 
Therefore, nine performance chart groups were reduced to analy- 
tical expressions and implemented on the TI-59 HHC. To show 
further possibilities and feasibility of implementation of the 
algorithms, they were 1) fully implemented on the Hewlett 
Packard-9830 (HP-9830) desk computer, 2) demonstrated on a 
microprocessor (INTEL Corporation Microcomputer System-48), 
and 3) considered for implementation on the A-7E onboard digi- 
tal computer and a microprocessor utilizing a recently developed 
number processing chip by the National Semiconductor Corporation 
(MM57109) . 
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II. DEVELOPMENT 



A. GUIDELINES 

The scope of this investigation was established after a . 
firm set of guidelines was defined. 

Being the official United States Navy standard for the 
A-7E aircraft, the A-7E NATOPS Manual was the sole source of 
performance data used to develop the algorithms. As such, 
and being subject to changes during the aircraft's life cycle, 
the need for possible future updates to the algorithms was 
acknowledged. The effective date of the NATOPS Manual from 
which these algorithms were developed is March 1975. Since the 
performance data yielded by the algorithms was identical to 
NATOPS Manual performance curves, the same restrictions and 
limitations apply. For example, takeoff airspeed calculation 
restricts the NATOPS Manual user to trailing edge flap positions 
between 20 and 40 degrees down (Figure 1). For that reason, one 
could not expect to calculate the flaps up takeoff airspeed 
using the developed algorithms . An additional feature provided 
by the algorithms was higher order interpolation. While the 
inexperienced NATOPS user might attempt to interpolate linearly 
between non-linearly spaced curves, the algorithms do not. 

An important guideline for the user's benefit was to ensure 
the execution of these algorithms after implementation was 
simple enough so very little training was required for the 
users. Intended users were Naval Flight Officers and Aviators. 
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CONDITIONS: 

MODEL: A-7E MILITARY RATED THRUST ENGINE; TF41.A-2 

DATA BASIS; FLIGHT TEST LANDING CONFIGURATION FUEL GRADE: JP-5 

DATE. NOVEMBER 1971 LEADING EDGE FLAPS DOWN FUEL DENSITY: 6.8 LB/GAL 



NAVAIR 01-45AAE-1 




SVOI - A1IOOT3A 



11-21 
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Takeoff Speed NATOPS Chart 



Not included in the scope of this thesis is an introduc- 
tion to the TI-59 HHC, HP-9830 desk computer and the INTEL 
Microcomputer System-48; however, to follow the computer pro- 
grams written for these devices would required their basic 
understanding . 

Another guideline established was that the performance 
calculators be light and small enough to be physically suited 
for its environment. For example, the TI-59 calculator and 
microprocessor could be used in a cockpit, briefing room or 
Air Operations Center. The HP-9830 desk computer would be 
restricted from cockpit use. 

Reliability was a necessary guideline. 

To make algorithm implementation on the TI-59 HHC feasible 
and since the program storing chip, the Continuous Read Only 
Memory (CROM), was limited to 5000 calculator program steps, 
the library of nine programs was required to fit into that 
space [3]. 

Finally, accuracy was a necessary consideration. The 
results obtained from the algorithms were required to be at 
least as accurate as following the performance charts manually. 
These accuracy requirements established were: One knox of air- 

speed, 100 feet of altitude or ground roll distance, 100 pounds 
of weight, ten seconds of time and one nautical mile of dis- 
tance . 
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B. PERFORMANCE CHART REDUCTION 

The reduction of the NATOPS Manual performance curves into 
analytical expressions was accomplished by a historically proven 
mathematical procedure, "least squares curve fitting". This ■ 
method was applied to certain A-7E performance data by LCDR 
W.M. Siegel (see Introduction, Section I). His brief explana- 
tion of the "Least Squares Fit Approximation (LSFA)" is included 
in Appendix A. 

Many performance charts from the NATOPS Manual contain 
three variables (two independent, one dependent) and are de- 
picted as a two-dimensional space with the third dimension 
illustrated by a family of curves. The reduction of such a 
chart can be accomplished as follows: 

1. Determine order of curves in family (i.e, second order, 
(y = A^ + A 2 X + A^x^). 

2. Apply LSFA to every member of the family of curves. 

3. Since the order of the curve families may vary, a 
general curve family could be depicted as follows : 



y = A, , + A, ^x + A, .x*^ + 



‘11 



12 ‘ 
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y = A^, + A^^x + A^oX^ + 
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( for curve z ) 
m 



‘ 21 ^ 



y=A, +A„x+A -x^ + . . 

^ ml m2 m3 

Apply LSFA to the coefficients. For example, plot 
. .A^^ versus z^, Z 2 » z^, respectively, yielding 



Doing the same with all coefficients 
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Given z and x, y can now be calculated by: 



Computing coefficients from equations generated in 



Step 4 . 

b. Applying coefficients to y = A^ + + . . .A^x^ 

6. It is important to note that although all curve family 
members must be of identical order, the equations representing 
the coefficients as a function of "z" need not be of similar 
order . 



Although applying LSFA to the family of curves and then 
to their coefficients was the normal method of chart reduction, 
it was not always used for the following reasons : 

a. Some charts were two-dimensional (LSFA still used). 

b. Some charts were reduced by inspection. 

(1) Linear curve families with linear spacing. 

(2) Time, distance, speed charts (d = v/t). 

c. Algorithm anomalies (see "Development Difficulties"). 

When used, the LSFA was accomplished by a program pre- 
written by the Hewlett Packard Corporation for use with the 
HP-9830. This program, although greatly facilitating the 
development portion of this investigation, was written for a 
two-dimensional problem and had to be executed at least once 

for each curve and once for each set of coefficients. 

A listing of all of the equations making up the perfor- 
mance algorithms are contained in Appendix C. The A-7E 
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performance chart groups from which they were developed are 
contained in Appendix B. They are in order: 

1. Low Level Cruise Performance. 

2. Takeoff Ground Roll Distance. 

3. Maximum Range Cruise Time and Speed at Constant 
Altitude . 

4. Maximum Range Cruise Fuel Required at Constant 
Altitude . 

5. Maximum Range Climb Airspeed Schedule. 

6. Takeoff Airspeed. 

7. Maximum Refusal Airspeed. 

8 . Optimum Endurance Altitude . 

9. Cruise Ceiling. 

Future reference in this thesis is made to algorithms and 
programs by the numbers above. 

C. EXAMPLE OF CHART REDUCTION 

An example of the procedure discussed in the previous sec- 
tion is presented below. The chart chosen for reduction is the 
lower graph of Figure 2, from Phase II of the A-7E Cruise Per- 
formance chart group . 

By inspection, all A^ and coefficients are equal to 
zero. The curves appear parabolic and therefore second order, 
yielding N = A^M^ . The example follows: 

N = intermediate result 

M = mach number 

D = drag count 
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CRUISE PERFORMANCE (A-7E) 



PHASE II - AIRCRAFT REFERENCE NUMBER 

MODEL: A-7E 

DATA BASIS; FLIGHT TEST 
DATE: NOVEMBER 1971 



n-117 



ENGINE; TF41-A-2 
FUEL GRADE: JP-5 
FUEL DENSITY: 6.B LB/GAL 






n-58 



Figure 2 



Cruise Performance 
Phase II NATOPS Chart 
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DRAG COUNT LINE 



CURVE EQUATION 



50 


N 


= 


1. 3915M 


10 0 


N 




2 . 7787M 


150 


N 


= 


4. 1658M 


200 


N 


z 


5. 5530M 


300 


N 


= 


8 . 3273M 


400 


N 


Z 


11.102M 



By plotting the A^ coefficients versus D (drag count), the 
LSFA yields: 

A^ = (4.3732E-3) + .027743D and therefore, 

N = ((4.3732E-3) + (. 027743D)M^ . 

This was a particularly simple chart to reduce but illus- 
trates the procedure. 

D. DEVELOPMENT DIFFICULTIES 

The normal method of reducing performance curves did not 
always yield useful information. One reason was although the 
NATOPS Manual Performance curves were constructed from experi- 
mental data, families of curves occasionally had very unusual 
spacing. They also were not always a true curve family; that 
is, they were of varying order. This can be visually detected 
in the lower graph of Phase III of the A-7E Cruise Performance 
chart group (Figure 3). The unequal and varying spacing be- 
tween curves with different "reference numbers" is obvious. 
Although the coefficients for each curve can be calculated, 
the coefficients determined for a LSFA equation for an inter- 
mediate curve would be incorrect. To be usable for the normal 
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CRUISE PERFORMANCE (A-7E) 



PHASE III - POUNDS OF FUEL PER NAUTICAL MILE 



ENGINE: TF4t*A'2 
FUEL GRADE; JP-5 
FUEL DENSITY; 6,8 LB/GAL 



MODEL: A-7E 

DATA BASIS; FLIGHT TEST 
DATE: NOVEMBER 1971 




Figure 3 



n-59 



Cruise Performance 
Phase III NATOPS Chart 
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method of chart reduction, a chart must have equal, constantly 
increasing, or constantly decreasing spacing between curves. 
When such an incompatible chart was encountered, it was neces- 
sary to interpolate between them. Two chart groups eliminated 
from consideration, "Fuel Required" and "Time to Climb from Sea 
Level to Selected Altitude", contained so many such curves 
(11), that very high order expressions would have been re- 
quired to compute the coefficients, making implementation on 
the TI-59 HHC impractical. The A-7E Cruise Performance lower 
chart of Phase I had the same anomaly (Figure 4). Because of 
the importance of the low level mission, however, the algorithm 
for this chart was developed, for sea level only though. The 
multiple algorithm was not developed but could have been for 
implementation on a desk computer. 

Another reason a straight application of LSFA was not 
always appropriate was the uniqueness of the upper graph of 
Phase I of the A-7E Cruise Performance chart group (Figure 4). 
This chart requires entry from the lower chart. A line is 
traced upward until the user contacts the appropriate Drag 
Count Line (dotted lines). The first pass through the Mach 
Number axis, a result of the lower chart, was defined M* . 
Instead of now tracing horizontally to the Transfer Scale axis 
(this value defined TS*), one must trace "between the solid 
guidelines" to the interception with a line traced vertically 
upward from the desired Mach number, M. The Transfer Scale 
would now be manually obtained by tracing horizontally to the 
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NAVAIR 01 45AAE-1 



CRUISE PERFORMANCE (A-7E) 



PHASE I - CLEAN AIRPLANE TRANSFER SCALE 

MODEL: A-7E ENGINE: TE4t-A-2 

DATA BASIS: FLIGHT TEST P'JEL GRADE: JP*5 





Figure 4 ii-s7 

Cruise Performance 
Phase I NATOPS Chart 
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TRANSFER SCALE 



vertical axis- To develop the algorithm for this problem, 
the equations of the guidelines were also calculated as a func- 
tion of Mach number. The values of the Transfer Scale resulting 
from M* intercepting the guidelines and tracing horizontally 
to the vertical axis were called TS^*, TS2*, ...TS^S from top 
to bottom. The original position, (M*, TS*), could now be 
determined in relation to (M*, TS *) and (M*, TS , ,*) • "n" and 

"n+1" indicate the upper and lower guidelines, respectively, 
which bracket (M“, TS*). This ratio provided the initial posi- 
tion relative to the guidelines: 



R = (TS*-TS *)/(TS*-TS ^ ’’O 
n+1 n n+1 

Using the desired Mach number, M, the Transfer Scales for 

the same two enclosing guidelines were calculated (TS^ and 

TS^_i_^) . The final position relative to the guidelines was 

maintained using the original ratio by solving: 

R = (x-TS ^,)/(TS -TS for x. 

n+1 n n+1 

"x" is the Transfer Scale with which the user now proceeds to 
Phase III of this performance chart group. Figure 5 depicts 



this problem graphically. 



E. ACCURACY 

A large number of results comparisons between the generated 
algorithms and manually traced performance problems were made. 
An infinite number of comparisons would be required to check 
all possibilities, but since the mathematical theory was so 
basic, the number of checks accomplished were considered suf- 
ficient . 
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Figure 5 

Guideline Chart Solution 



TRANSFER SCALE 



All nine algorithms were checked for accuracy on the 
HP-9830 desk computer. The number of checks for each algo- 
rithm was proportional to the ease of manually tracing through 
the performance charts. The author spent considerable time 
obtaining performance results from the NATOPS Manual charts 
and a relatively small amount of time computing the problems 
on the desk computer once the algorithms had been implemented. 

In a significant number of instances, the results disagreed, 
but after rechecking, the solution obtained manually was in 
error. This supported the contention that manual manipulation 
of the performance charts is an error prone procedure, even 
with an experienced user. 

In a few rare instances , the author entered the required 
given data incorrectly into the desk computer. These miskeying 
errors , not procedural , were noticed as soon as the answer was 
produced. A user familiar with the A-7E performance character- 
istics would normally notice an answer resulting from grossly 
incorrect data input. It is acknowledged, however, that there 
is no failsafe check on the programs. When using a desk com- 
puter, the required input data can be printed along with the 
answer to ensure the user of the correctness of the input data. 
For a hand held calculator, however, computing a performance 
problem twice would provide a check, which is what many NATOPS 
Manual users often do. As with all computer programs, a desired 
result requires accurate input data. 
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Except for those noted below, the results of programs 
checked (using five significant figures) were indistinguish- 
able from the answers obtained by manually manipulating the 
performance charts . Answers produced from the algorithms 
were rounded off to the nearest digit. 

PROGRAM MAXIMUM DEVIATION 

Maximum Refusal Speed 2 knots 

Takeoff Airspeed 1 knot 
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III. IMPLEMENTATION 



A. DESK COMPUTER 

The use of a desk computer capable of producing A-7E per- 
formance information within seconds (less than three seconds 
computation time for the longest algorithm) would be ideal for 
a squadron briefing room or Air Operations Center use. The 
HP-9830 desk computer was used for this implementation stage. 

Very little training would be required for personnel to load 
the programs stored on a cassette tape cartridge and execute 
them. 

A knowledge of "basic" computer language is required to 
fully understand the nine HP-9830 programs in Appendix D [4]. 

The nine programs are in the same order as the algorithms of 
Appendix C. 

Only in the Low Level Cruise Performance program are sub- 
routines required for linear interpolation or for the iterative 
method to find the Transfer Scale (see "Development Difficulties") . 
All other programs are straight forward, sequential computations. 
In these programs, the coefficients defining a curve (y = f(x)) 
for a given set of conditions are calculated. That chart 
result, "y", is then calculated for the given independent vari- 
able "x". The next chart of that group is similarly treated 
and so on until the "final result" is achieved. 

The HP-9830 programs are very useful since they prompt the 
user to supply the correct information. Most of the programs 
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"request, then accept" those inputs required for the applicable 
NATOPS Manual performance chart. The HP-9830 then prints the 
data just entered (ensuring the user that data input was as 
desired) followed quickly by the solution. The computer is 
instantly ready to receive new data for another calculation. 

Programs 1, 2, 3, 4 and 7 (as identified in "Performance 
Chart Reduction, section II-B), are written in this "request, 
then accept" format. The shorter programs, 5, 6, 8 and 9, 
were written with an initial set of input data already in the 
program. This format allowed the computer to step incrementally 
through the allowable range of values for the input data, thus 
calculating a "table of performance data" for the applicable 
performance chart group. These programs are easily altered to 
the "request, then accept" format by some simple edit commands 
[4]. 

The variables used in the programs are defined following 
each program in Appendix D. 

B. HAND HELD CALCULATOR 

The many favorable features of the hand held calculator 
encouraged its implementation of the performance algorithms . 

Its small size allowed consideration for use in the cockpit. 

Its simplicity and reliability was an advantage making it 
especially suited for users of varying experience (including 
no experience). Although its execution speed was the slowest 
of all devices used, the computation time was still much faster 
than using the NATOPS Manual. 
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The Texas Instruments- 5 9 (TI-59) programmable hand held 
calculator (HHC) was selected for implementation. This selec- 
tion was made for several reasons. At the time, it was the 
only calculator available to the author which allowed permanent 
program storage (on magnetic cards). Additionally, the Texas 
Instruments Corporation had the capability to combine all pre- 
written performance programs, up to a 5000 program step limit, 
onto a Continuous Read Only Memory (CROM) chip, making the 
A-7E performance programs a permanent part of the calculator. 
This CROM chip can also be used on the less expensive TI-58 HHC. 
These features made the TI-58/59 (with CROM) a practical system 
for the A-7E Naval Aviation community. 

One might consider the calculator's inability to prompt the 
user for inputs a shortcoming of this implementation candidate, 
but a company spokesman, Mr. Richard Cuthbert, stated a new 
face could be fitted onto the calculator, identifying different 
buttons with the input data categories such as GW for gross 
weight, FLPS for flap position, T(°C) for temperature, and so 
on C 3 ] . 

Some time was required for the author familiarization with 
the TI-59 HHC and its capabilities. For a detailed explanation 
of comments in this section involving TI-59 programming and 
Appendix E, consult the user's manual [5]. 

All programs were entered with the calculator memory parti- 
tioned to allow 879 program steps and ten memory storage 
locations. The loss of program steps in order to provide 
coefficient storage locations (ten to one) was the reason for 
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partitioning in this manner. Only five significant figures 
were considered necessary for computational accuracy. Consider- 
ing the number possibilities (1.2345 to 1.2345E-12) might take 
from six to ten program steps, this was less than the absolute 
ten program steps sacrificed for a storage location. The ten 
memory storage locations were used to store the input data at 
program execution start but were often reused after the input 
data storage was no longer required. 

The programming language level of the TI-59 HHC is below 
the HP-9830 's and above a microprocessor's (discussed later) 
in sophistication. The algorithms were computed in a more 
space-saving manner than on the HP-9830. For example, in com- 
puting a first order polynomial, the HP-9830 program functioned 
as follows: 

B(0) = 

B(l) = + A 22 Z 

y = B(0) + B(l)x. 

The TI-59 HHC was programmed to compute as follows: 

(All “^ 12 ^^ ^"^21 ^ 22 ^^^ ■ 

In the Low Level Cruise Performance program, the linear 
interpolation and iterative methods to follow guidelines (dis- 
cussed in previous section) was still accomplished using the 
more tedious TI-59 HHC language. 

Using the partitioning already described, a program limit 
of 879 program steps was imposed (filling two magnetic cards). 
Two programs, "Takeoff Ground Roll Distance" and "Low Level 
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Cruise Performance", exceeded this limit and had to be continued 
on extra cards. These programs were written to allow storage 
of an intermediate result into the T-register. The rest of the 
cards could then be read in, any lost or newly acquired input 
data entered, and program execution would continue, automatically 
retrieving the stored intermediate result from the T-register. 
These artificial necessities for program completion using the 
magnetic cards would not be necessary if the programs were 
stored permanently in the CROM. 

The total number of steps required for the nine performance 
algorithms programmed on the TI-59 HHC was 5461 steps. By sub- 
routining (340 steps of programming are common to two programs), 
the total number could be reduced to 5121 steps. The elimina- 
tion of the artificial steps required for the oversized programs 
would reduce the overage more. The sole intent of this imple- 
mentation phase was not to fit these nine programs into the 
5000 step CROM. If the inclusion of all nine programs was 
desired, streamlining aid offered by engineers from the Texas 
Instruments Corporation plus the reduction of significant 
figures in a non-critical area would accomplish this. 

The program listings, storage location usage, user instruc- 
tions, and execution times are included in Appendix E. 

C. MICROPROCESSOR 

1 . Single Board Computer using Software for Mathematical 
Operations 

The single board computer (SBC) implementation was 
investigated both as an extension of thesis work and to meet 
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the course objectives of AE-4900, Air Data Systems. Work 
toward this effort was also done by LCDR W.M. Siegel. The 
performance algorithms were to be processed on a SBC using an 
INTEL Corporation 8048 Programmable Read Only Memory (PROM), 
external random access memory (RAM) and a program counter. 
Software development was completed on the INTEL Prompt-48 
(Microcomputer System-48 language) using an INTEL 8035 arith- 
metic logic unit (ALU) . Although a SBC using the 8048 PROM 
and requiring a digital keyboard and display was never actually 
constructed because of the time limitations , the software opera- 
tion was successfully demonstrated on the Prompt-48. 

To preserve the programs between operation periods, 
the Prompt 48 was hand wired as specified in the user's manual 
to an ASR-35 Teletype set which allowed paper tape storage [6]. 
The Prompt— 4 8 provided 10 2 4 by two bytes of RAM and 6'4 by two 
bytes of resident memory. Although the MCS-48 instruction set 
will not be discussed in this thesis, a basic understanding of 
assembly level language is necessary to understand the developed 
software presented in Appendix F [7]. This microprocessor 
program listing includes the MCS-48 instructions in hex code 
and literal mneumonics and includes full documentation to facili- 
tate interpretation. 

A full performance algorithm was not implemented on the 
Prompt-48 because of its memory storage limitations. The 
original intent was to exercise the software of the complete 
A-7E Takeoff Ground Roll Distance algorithm on the Prompt-48. 
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After the necessary routines were written and stored, only room 
for three coefficients remained (98 coefficients required for 
this algorithm) . Since implementation capability was the 
desired result, the computation of a second order polynomial 
was considered sufficient. Although this effort was software 
oriented, the necessary RAM storage for the additional coeffici- 
ents and executive routine could have been easily provided for 
a SBC, 

The software development for algorithm implementation 
required routines for input/output (I/O), executive direction, 
binary to binary coded decimal (BCD) and BCD to binary conver- 
sions , and floating point binary addition and multiplication 
routines. The I/O and executive routines were written by LCDR 
Siegel. The nonavailability of a number oriented microprocessor 
at the time of this effort required the development of the 
mathematical package described above. The advantages for such 
a capability will be discussed in the following section. 

In addition to the microprocessor software developed by 
the author and LCDR Siegel, the I/O and display routines would 
require alteration for SBC implementation since a digital dis- 
play and keyboard would replace the Prompt-48. 

Figure 6 illustrates the solution method. Figure 7 is 
a flow chart of the program execution sequence. Figures 8 and 9 
show the Prompt-48 RAM and resident register memory, respec- 
tively. 
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SBC SOLUTION METHOD 




Figure 6 

SBC Solution Method 
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SBC PROGRAM EXECUTION SEQUENCE 




Figure 7 

SBC Program Execution Sequence 
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RANDOM ACCESS MEMORY MAP 



ADDRESS 
000-069 
06 A- 06 F 
070-079 
07A-0C6 
0C8-0E2 
0E5-0FF 
100-2EC 
300-3FF 



USE 

INPUT AND DISPLAY 
EXECUTIVE ROUTINE SEGMENT 
COEFFICIENT STORAGE 
MAIN EXECUTIVE ROUTINE 
BINARY TO BCD EXECUTIVE ROUTINE 
MISCELLANEOUS SUBROUTINES 
ADDITION AND MULTIPLICATION SUBROUTINES 

BCD TO BINARY EXECUTIVE ROUTINE AND CONVERSION SUBROUTINES 



Figure 8 

Random Access Memory Map 
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RESIDENT REGISTER MAP 



ADDREE 


;s USE 


ADDRESS 


USE 


20 


LSB 


30 




21 


X- LOCATION 


31 




22 


ARITHMETIC REGISTER 


32 


LSB 


23 


MSB 


33 


DISPLAY HEX 


24 


EXPONENT 


34 


BUFFER 

j 


25 


LSB Y- LOCATION | 


35 


MSB j 

i 


26 


MSB ARITHMETIC REGISTER 


36 


DECIMAL POINT MASK j 

! 

. __ i 


27 


EXPONENT 


37 


CHARACTER COUNTER j 

. . ... . j 


28 


LSB BCD- BINARY 


38 


LSB 1 


j 29 

r 


MS B CONVERS ION 
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I 

j 

i 

t 


i 

2A 


EXPONENT 


3A 


DISPLAY j 

1 


2B 




3B 


1 


2C 




3C 


BIT 


2D 




3D 




2E 




3E 


PATTERNS 


2F 




3F 


MSB 



Figure 9 

Resident Register Map 
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The second order polynomial, y = B(0) + B(l)x + B(2)x^, 
was calculated using a mathematical executive routine (alter- 
able for any size polynomial and any number of polynomials). 

The only mathematical operations required were multiplication 
and addition of positive or negative numbers. For speed, binary 
arithmetic was used. For increased storage capability and 
mathematical efficiency, a floating point capability was in- 
cluded . 

The calculation routine proceeded as follows : 

B2*x = (B 2 X) 

(B2X) + (Bj_) = (B2X + B^) 

(B 2 X + Bj^)*x = (B2X^ + B^x) 

(B2X^ + B^x) + Bq = (B^x^ + + Bq) 

Although all mathematical operations are performed in 
the 8-bit (2-byte) accumulator register of the 8035 ALU (for 
a SBC, the 8048 PROM), a working accumulator using five regi- 
sters (resident memory registers two through six) , was 
established. All numbers in the program (independent variable 
"x" after conversion to binary, coefficients stored in RAM 
070-079 and the 'result’) were in one of two binary conventions. 
While in storage, the numbers were in "storage" convention. 

The numbers were shifted from "storage" to "working" convention 
only when transferred from the X and Y locations (see resident 
register memory map. Figure 9) to the working accumulator 
(registers two through six) . When the desired operation was 
completed, the result was returned to the "storage" convention 
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and moved to the "X" location. Figure 10 displays the "storage 
and working" conventions . 

This software was successfully demonstrated on the 
Prompt-48. The user instructions for the Prompt-48 to repeat 
the demonstration are listed below: 

(1) Ensure the 8035 ALU or 8048 PROM is inserted in the 
"execution" socket of the Prompt-48. 

(2) Enter the program in hex code in the proper storage 
locations as listed in Appendix F. 

On the Prompt-48, press the following keys to clear 
the resident register memory: 

"C" 

"Registers" 

II 0" 

It It 
) 

11411 

" 8 " 

Do not press "Program Memory" instead of "registers" 
or the program just entered will be erased. 

(3) To execute the program, press the following keys: 

"A" 

IT 2 

"Execute" 

"Go" 

"No Break" 

" 0 " 

"Execute" 
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BINARY CONVENTIONS 



Each large block depicts 1 byte which includes 4 bits. The compart- 
mented blocks represent 1 bit. 



STORAGE CONVENTION 





t i S 


i 


\ 1 


L_ 


* ■ 



EXPONENT 



t I r 6 



Jl. 



r 

L_L 



■7~r 



-MANTISSA. 



*o 1 / 



T— 

1 


1 

L 1 


L« 


— I L_ 



• ^ is 

T T 

I 



"T 



DECIMAL 
SIGN BIT 



X 



TZ o (g It I 
— < — T — H- 

> ' I 



In storage convention, the mantissa is left justified to bit 5, A 
positive number is denoted by 0 in the first bit of the second byte (sign 
bit); a 1 indicates a negative number. 



WORKING CONVENTION 




In working convention, the mantissa is left justified to bit 4. 

The sign bit is stored in FO (X-location number) and FI (Y-location num- 
ber) flags of the program status word. 



Figure 10 

Binary Conventions 
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The display will blank, awaiting the input of the 
independent variable "x" . To enter "x" , enter the digit keys 
for numbers (base 10) and "D" for decimal point. "x" will be 
displayed on the digital display as it is entered. To compute 
the algorithm (second order polynomial), press "E". The answer 
will rapidly appear. To calculate the polynomial with a new 
value for "x", start at Step 3. 

(4) To prevent the time consuming reloading of the pro- 
gram, it is advisable to store the program on a peripheral 
device (paper tape, disc, etc.). 

2 . Single Board Computer using Number Oriented Micropro- 
cessor 

Very recently, the National Semiconductor Corporation 
began production of a chip intended for use in number processing 
applications [8]. This chip, the MM57109 MOS/LSI, is capable 
of all scientific calculator functions, test and branch capa- 
bilities, internal number storage, and I/O instructions. Of 
the specific calculator functions, only addition, subtraction 
and multiplication would be used. 

A SBC using this chip would need the 8048 PROM for 
coefficient and executive routine storage but would not need 
the space consuming mathematical package of the SBC in the last 
section. A program counter would still be required but external 
RAM would not. The computation time would be increased over the 
demonstrated SBC (approximate computation time of a HHC), but 
the simplicity of programming would make this proposed SBC very 
attractive . 
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D. A-7E TACTICAL COMPUTER 



In February 1978 the author made a trip to the Naval Air 
Facility at China Lake, California. The purpose of this visit 
was to receive indoctrination on the TC-2/2A tactical computers 
and obtain a programming manual for these devices. The desired 
goal was implementation of selected performance algorithms on 
the laboratory bench computer run by the A- 7 Program Office of 
the Naval Weapons Center (NWC). A thorough understanding of 
the computer's capabilities and limitations was provided by Mr. 
Robert Westbrook, a software technician. 

The A-7E computer provides very accurate navigation and 
weapons guidance capability. The TC-2 and TC-2A computers are 
a generation apart, the TC-2A being over two times faster and 
having twice the storage capability of its earlier version. 

Both computers are operational at this time. Specific design 
and programming information is available from the programming 
manual [ 9 ] . 

The instruction set of the tactical computer provides fixed 
point arithmetic, logical transfer of control (branching), 
address modification and single word input/output instructions 
specifically intended for operations primarily involving arith- 
metic. These features made the implementation of algorithms a 
logical decision. Several factors made this implementation by 
the author impractical. The computer design was quite old, the 
instruction set being very tedious and difficult to interpret. 
The computer's inability to function using floating point 
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arithmetic would require a significant software effort in that 
area alone. The time required to become fully familiar with 
the instruction set, write the software, and load and test the 
programs at NWC would have been prohibitive for this investiga- 
tion. 

It is hoped that the programmers at NWC will be able to 
implement those algorithms deemed desireable to achieve an 
onboard capability. Takeoff Airspeed and Maximum Refusal Air- 
speed are considered ideal for implementation. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



Nine of the A-7E NATOPS Manual performance chart groups 
were reduced to a series of analytical expressions or algo- 
rithms. These algorithms, accurate to five significant figures, 
are as accurate as results obtained by manual manipulation of 
the performance charts . 

Implementation was made on three data processors of different 
programming levels and storage capabilities. These devices and 
degrees of implementation were: 

(1) HP-9830 Desk Computer - complete implementation 
with successful demonstration. 

(2) TI-59 Hand Held Calculator - complete implementation 
with successful demonstration. 

(3) Microprocessor - partial implementation with suc- 
cessful demonstration. 

In view of the success of this investigation, recommendations 
concerning implementation possibilities are listed below: 

(1) Complete reduction of the NATOPS Manual performance 
charts could be accomplished and implemented onto a desk com- 
puter as one large program capable of performance data computation 
within seconds. The desk computer would be ideal for mission 
planning on a squadron or air wing level or for Air Operations 
Center use. 

(2) The programs written for the TI-59 HHC could be con- 
solidated onto a CROM and used with a TI-58 HHC for use on a 
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squadron level. As an alternative, the software could be 
rewritten for any HHC of comparable capability. 

(3) Although implementation on a single board computer 
using a number oriented microprocessor is completely feasible, 
because of programming ease and cost consideration, the HHC is 
considered a superior implementation possibility at this time. 

(4) The A-7E tactical computer could easily be pro- 
grammed by software engineers at NWC, China Lake, California, 
to produce an onboard capability. 
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APPENDIX A 



Least Squares Fit Approximation 

References 10 and 11 describe the Least Squares Fit Approxi- 
mation in detail. In general the problem is to represent a set 
of "n" data points in two-dimensional space 

X^, i = 1 to n 

by a polynomial expression of a curve whose degree is less than 
"n" . Two classes of problems exist: 

(1) Linearly independent - those in which the degree 
(d) of the polynomial is one less than the number of data 
points 

d = n-1 (1) 

(2) Linearly dependent - those in which the degree (d) 
is less than n-1 

d < n-1 (2) 

As an example, a set of four (4) data points randomly spaced 
was chosen. If a third degree polynomial of the form 

Y = A + BX + CX^ + DX^ (3) 

were desired, and the data points X^ and Y^ were inserted 
(i = 1 to 4) into four such equations, an exact solution for the 
four unknown coefficients would exist . These four unknowns 
could be found from the four equations by numerous conventional 
techniques (Direct substitution, Cramer's rule, etc.). The poly- 
nomial expression generated would be termed a "col-location" 
polynomial because its plot would pass through all data points. 
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It is often advantageous to describe a set of data points 
by a curve that does not pass through each point. This type 
of polynomial would be termed a "regression" equation. For any 
set of data points an infinite number of regression expressions 
exist for any specified degree (except the linearly independent 
case) and the object of the Least Squares Method is to find the 
polynomial coefficients of the chosen degree that best describe 
the data points. In the previous example of four data points, 
assume that, instead of the third degree form chosen, a second 
degree equation were selected of the form 

Y = A + BX + CX^ (4) 

With four data points, the polynomial is overspecified and thus 
linearly dependent. For this case an infinite number of solu- 
tions exist for the coefficients a, b and c. If an error term 
(5) were defined for any given X,Y pair as 

6i = lY^ - A + BX^ + CX^^I (5) 

a total squared error term (E) could then be defined by squaring 
and Slimming the terms attained: 

N 

E = 1 6 2 (6) 

i = l 

If E were them minimized for any given degree chosen, the best 
Least Squares Fit would have been achieved. 

If the values for 6 from Equation 5 were inserted in Equation 
6 and the partial derivative of E were taken with respect to the 
coefficient A, an equation would be generated that when set equal 
to zero (0) would define a minimum value of E for a given value 
of A. If the same operation were performed with respect to the 
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coefficients B and C then three equations would be generated 
with three unknowns (A, B and C) . The solution of these simul- 
taneous equations would produce the coefficients A, B and C, 
that would minimize the value of E and hence would produce a 
Least Squares Fit approximation to a set of linearly dependent 
equations . 

A numerical procedure has been developed to accomplish this 
task. An example of this procedure has been included in the 
following paragraphs [10, 11]. 

Least Squares Fit Method Example 
Given the following set of data: 

X 0 1 2 4 7 

f(X) = Y 0 1 3 12 20 

fit a curve of the form 

f(X) = Y - A + BX + CX^ 

STEP 1: Substitute all pairs of data into the form equation 

yielding the fact that the coefficients (A, B and C) must 
satisfy all the following: 

0 = A + B(0) + C(0)^ 

1 = A + B(l) + C(l)^ 

3 = A + B(2) + C(2)^ 

12 = A + B(4) + C(4)2 

20 = A + B(7) + C(7)^ 

Now multiply each expression by its coefficient of A in that 
expression and add all equation yielding 
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36 



5A + 14B + 70C 



Now multiply each expression by its coefficient of B in that 
expression and add all the equations yielding 



0 = 


OCA) 


+ 0(B) 


+ 0(0 


1 = 


A 


+ IB 


+ 1C 


6 = 


2A 


+ 4(B) 


+ 8(0 


48 = 


4A 


+16(B) 


+ 64(0 


140 = 


7A 


+44(B) 


+ 343(0 


195 = 


14A 


+70(B) 


+ 416(0 



Now multiply each expression by its coefficient C in that 
expression and add all the expressions yielding 

0 = OCA) + 0(B) + 0(C) 

1 = 1(A) + KB) + 1(C) 

12 = 4(A) + 8(B) +16(0 

192 =16(A) +64(B) +256(0 

980 =49(A) +343(B)+2401(C) 

1185 = 70A + 416B + 2674C 

Now solve the following three previously generated equations 
the coefficients A, B and C yielding 
36 = 5A + 14B + 70C 
195 =14A + 70B +416C 
1185 =70A +416B +2674C 

A = -.99, B = 2.6, C = .065 

and 

Y = -.99 + 2.6X + .065X^ 



for 
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The following plot and chart depict the original data and 



the data 


obtained from the 


equation for the fitted 




Original 


Fitted Curve 
Polynomial 


X 


Y 


Y 


0 


0 


-.98 


1 


1 


' 1.67 


2 


■5 


4.48 


4 


12 


10.46 


7 


20 


20.41 



Q . E . D . 
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APPENDIX B 



NATOPS Manual Performance Charts 

These charts from which the performance algorithms were 
developed are listed below in order: 



Figure 


Title 


B1 


Cruise Performance, Phase I 


B2 


Cruise Performance, Phase II 


B3 


Cruise Performance, Phase III 


B4 


Cruise Performance, Phase IV 


B5 


Takeoff Factor 


B6 


Takeoff Ground Roll Distance 


B7 


Adjusted Takeoff Ground Roll Distance 


B8 


Maximum Range Cruise at Constant Altitude (Time, Speed) 


B9 


Maximum Range Cruise at Constant Altitude (Fuel 
Required) 


BIO 


Military Power Climb Schedule 


Bll 


Takeoff Speed 


B12 


Maximum Refusal Speed 


B13 


Cruise Ceiling and Optimum Endurance Altitude 



50 



NAVAin 01-4GAAE-1 



CRUISE PEREORMAtlCE (A-7E) 



PHASE I - CLEAN AIRPLANE TRANSFER SCALE 

MOOeU: A-76 ENGINE: TF4I.A-2 

DATA BASIS: FLIGHT TEST PUCL GRAOE: JP-5 

DATE: NOVEMHUR 1971 FUEL DENSITY: 6.0 LO/GAL 





Figure B1 

Cruise Performance, Phase I 
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NAVAIR 01-45AAE1 



CRUISE PERFORMANCE (A-7E) 



PHASE II - AIRCRAFT REFERENCE NUMBER 

MODEL: A-7E 

DATA BASIS: FLIGHT TEST 

DATE: NOVEMBER 1971 



11-117 



ENGINE: TF41-A.2 
FUEL GRADE: JP*5 
FUEL DENSITY: G.8 LO/CAL 




Figure B2 

Cruise Performance, Phase II 
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NAVAIR 01-45AAE1 



CRUISE PERFORMANCE (A-7E) 



PHASE III - POUNDS OF FUEL PER NAUTICAL MILE 



MODEL: A-7E ENGINE: TF41.A-2 

DATA OASIS: FLIGHT TEST FUEL GRADE: JP-5 

DATE: NOVEMoen 1971 FUEL DENSITY; 6.8 LQ/GAL 






Figure B3 

Cruise Performance, Phase III 



n-59 
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NAVAIR 01-45AAE-1 



CRUISE PERFORMANCE (A-7E) 



PHASE IV - FUEL FLOW 

ENGINE: TF41-A-2 
FUEL GRADE: JP-5 
FUEL DENSITY: 6.B LB/GAL 





MODEL: A-7E 

DATA BASIS; FLIGHT TEST 
DATE: NOVEMBER 1971 



- 03 - 



11-60 



Figure B4 

Cruise Performance, Phase IV 
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NAVAIR 01*45AAE*1 



TAKEOFF FACTOR (A-7E) 



MOOEL: A-76 ENGINE: TF41-A-2 

DATA BASIS: FLIGHTIEST FUE L GRADE : JP-5 

DATE: NOVEMBER 1971 FUE L DENSITY: 6.8 LB/GAL. 







n-18 Changes 

Figure B5 

< 

Takeoff Factor 
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NAVAIR 01-45AAE-1 



I j TMEOFF GROUND ROLL DISTANCE (A-7E) 



MODEL; A-7E 

DATA BASIS: FLIGHT TEST 
DATE: NOVEMBER 1971 



CONDITIONS: 

LEVEL HARD SURFACE RUNWAY 
MILITARY RATED THRUST 
LANDING CDNFIGURATION 
ZERO HEADWIND 
CG: 26% MAC 
FULL FLAPS 



ENGINE: TF41-A.2 
FUEL GRADE: JP-5 
FUEL DENSITY: 6^ LB/GAU 



^NOTE^ 



For minimum ground roll corresponding to minimum 
liftoff speed, subtract 500 feet. 

For humidity effects on takeoff distance, ground roll 
distecKes should be increased 1% for each 10% increase 
In the relative humidity above 40%, 




76E287(l)-04- 74 



Change 6 11-19' 

Figure B6 

Takeoff Ground Roll Distance 
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NAVAIR 01-45AAE-1 



TAKEOFF GROUND ROLL DISTANCE (A-7E) 



ADJUSTED GROUND ROLL DISTANCE 

model: A-7E CONDITIONS: 

DATA BASIS: FLIGHT TEST HARD SURFACE RUNWAY 

DATE: NOVEMBER 1971 MILITARY RATED THRUST 

LANDING CONFIGURATION 
LEADING EDGE FLAPS DOWN 

For humidity affects on takeoff distance, ground roll 
distances should ba increased 1% for each 10% increase 
in the relative humidity above 40%. 



GROUND ROLL DISTANCE - 1.000 FT 





ENGINE; TF41-A-2 
FUEL GRADE; JP-5 
FUEL DENSITY: 6.B LB/GAU 



11-20 



Figure B7 

Adjusted Takeoff Ground Roll Distance 
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NAVAIR 01-45AAE-1 
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NAVAIR 01-45AAE1 



MAXIMUM RANGE CRUISE AT CONSTANT ALTITUDE (A-7E) 



FUEL REQUIRED 

MODEL; A-7E 

DATA BASIS: FLIGHT TEST 

DATE: NOVEMBER 1971 



ENGINE; TF41A-2 
FUEL GRADE; JP-5 
FUEL DENSITY: 6.8 LB/GAL 





FUEL FLOW 1.000 LB/HR 
0 2 : 4 6 8 10 





76E270(7)-03- 72 



Figure B9 
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Maximum Range 



Cruise at Constant Altitude (Fuel 



Required) 
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NAVAIR 01-45AAE-1 



MILITARY POWER CLIMB (A-7E) 



CLIMB SPEED SCHEDULE 

ENGINE: TF41-A.2 
FUEL GRADE: JP*5 
FUEL DENSITY: 6.8 LB/GAL 



MODEL: A-7E 

DATA BASIS; FLIGHT TEST 
DATE: NOVEMBER 1971 




76r2r6{ 72 



n-48 



Figure BIO 

Military Power Climb Schedule 
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NAVAIR 01-4GAAE-1 
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NAVAIR 01-45AAE-1 
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Maximum Refusal Speed 



NAVAIR 01-45AAE-1 



CRUISE CEULUG AHD OPTIMUM EHDURAIICE ALTITUDE (A-TE) 



MODEL; A-7E 

DATA BASIS: FLIGHT TEST 
DATE: NOVEMBER 1971 



ENGINE: TF41-A.2 
FUEL GRADE; JP-5 
FUEL DENSITY; 6^ LB/GAL 




15 20 25 30 35 40 45 

GROSS WEIGHT AT ALTITUDE - 1,000 LBS 




Figure B13 

Cruise Ceiling and Optimuin Endurance Altitude 



11-53 
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APPENDIX C 



Generated Algorithms 
LOW LEVEL CRUISE PROGRAM 



Phase I 

Ml = -92.512 + 236.896G 
Transfer Scale Versus Drag Count 

AO = -2.3287 -. 26316D + .0073327D^ - ( 7 . 513E-5 ) D^ + ( 3 . 5 39 6E- 7 ) D'* 
-(7. 78E-10)D^ + (6.462E-13)D® 

A1 = 4.835 + 1.0956D -.030653D^ + ( 3 . 1912E-4 )D^ - ( 1 . 5 2 7 6E-6 ) D^ 

+( 3 .408E-9)D® -(2.8692E-12)D® 

A2 = 10.284 -1.0719D + .031094D^ - ( 3 . 28 78E-4 ) D^ + ( 1 . 59 5E-6 ) D*^ 
-(3.6009E-9)D^ + ( 3 . 0 6 34E-12 ) D® 

SI = AO + (A1)(M1) + (A2)(M1)^ 

Transfer Scale Versus Guidelines 

BO = 22.819 -31.7341 + 41.331^ -5.09531® 

B1 =-154.98 + 217.511 -261.731^ + 35.9051® 

B2 = 405.08 -525.561 + 607.491^ -88.7371® 

B3 = -445.62 + 542.981 -611.551^ + 92.8941® 

B4 = 184.78 -204.421 + 225.891^ -35.1891® 

S = BO + (B1)(M1) + (B2)(M1)^ + (B3)(MD® + (B4)(M1)'^ 

Phase II 

R = S + 2C(4. 3732E-3) + .027743D]M^ 
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Phase III 



BO = 5.6253 -1.989R + 3.0252R^ -1.0761R^ + • 17675R'' -.013095R® 
+ (3.526E-4)R® 

B1 = 205.3012 -248.9317R + 91.66355R^ -15.55218R^ + 1.224432R'' 
-.0395333R® + ( 2 . 896385E-4 )R® 

B2 = -1052.123 + 1231. 24R -487.4233R^ + 91.6522R^ -8.662962R‘* 

+ .3953974R^ - . 006905535R® 

B3 = 1680.142 -1950. 139R + 788.8513R^ -152.5733R^ + 15.03819R'* 
-.7274139R® + .013707R® 



R3 = R 








R1 = 2 


(Integer (R/2) ) 




R2 = R1 


+ 2 






N1 = BO 


+ (B1)(R1) 


+ (B2)(R1)^ 


+ (B3)(R1) 


N2 = BO 


+ (B1)(R2) 


+ (B2)(R2)^ 


+ (B3)(R2) 



Using Linear Interpolation 
N = N1 + [(N2-N1) (R3-RD/2] 

P = 4.9746N + (7.9043E-6)N^ 

Phase IV 

N4 = [6.4375 + .010426T - ( 6 . 89 2 5E-6 )T^ + ( 4 . 912 7E- 7 ) ]M 
F = .1(N4)P 



TAKEOFF DISTANCE PROGRAM 

BO = 13.086 -.00017113A - ( 2 . 0 6 55E- 7 )A^ + ( 3 . 6861E-11 )A^ 
-(2.4156E-15)A‘' 



B1 = -.045635 -(7 . 8931E-6)A + ( 3 . 7545E-9 )A^ - ( 9 . 70 8 8E-1 3 ) A^ 

+ (6 .997E-17)A** 

B2 = -.001317 -(8.2558E-7)A + ( 4 . 0 7 39E-10 ) A^ -( 8 . 548E-14) A^ 

+ (5 .4964E-18)A‘* 

B3 = -(1.9097E-5) + (1.3671E-8)A - ( 9 . 469 4E-12 ) A^ + ( 2 . 0434E-15 )A^ 
-d.4617E-19)A‘* 

C = BO + Bl(B) + B2(B)^ + B3(B)^ 

If double datum on, 

E = 1.9773 + .56598C 
If double datum off, 

E = .54178 + .65876C 

GO = -(4.8896E + 5) + (8.4974E + DG - ( 5 . 785 6E- 3 ) G ^ + ( 1 . 9 37 3E-7 ) G ^ 

-( 3.1744E-12)G‘* + ( 2 . 0446E-17 ) G^ 

G1 = (5.8621E+4) - ( 1 . 0146E+1 ) G + ( 6 . 8 8 0 7E-4 ) G^ - ( 2 . 29 2E-8 ) G ^ 

+ ( 3.7387E-13)G‘* - ( 2 . 39 64E-18 )G^ 

H = GO + GKE) 

If relative humidity <40%, K = H 
If not, K = 4{[(I-40)/1000]+l} 

LO = 67.124 + .89509K + ( 2 . 3306E-5) - ( 1 . 6 2 54E-9 ) 

+ ( 3. 3728E-14)K‘* 

LI = -9.0995 -d.0856E-2)K + ( 2 . 1754E-7 ) - ( 2 . 5 3 2 7E-11 ) 

+ (1.197E-15)K‘* 

L2 = (1.4782E-1) - ( 2 . 166 6E-6 ) K + ( 3 . 4274E-9 ) - ( 2 . 7817E-13 ) 

+ (9 . 3077E-18)K‘* 

M = LO + LKL) + L2(L)^ 

If winds calm, M = K 
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XO = (4.5704E+1) + .93429M + ( 2 . 2 265E- 5 )M^ - ( 2 . 3 3 8E-9 ) 

+ (7.941E-14)M‘* 

XI = 7.9472 + .014914M + ( 9 . 0 708E-6 )M^ - ( 7 . 12 35E-10 )M^ 

+ (3 .0684E-14)M‘^ 

X2 = 5.3616 -.0085136M + ( 3 . 5914E-6 )M^ - ( 4 . 59 3 2E-10 ) 

+ d.9889E-14)M'* 

X = XO + XKN) + X2(N)^ 

QO = 2604.2 -2.1694X + .0010915X^ - ( 1 . 1119E-7 ) X^ + ( 3 . 6 6 2E-12 ) X‘* 
Q1 = -175.73 + .22601X - ( 7 . 5225E-5 ) X^ + ( 7 . 7018E-9 ) X^ 
-C2.5437E-13)X'* 

Q2 = 2.8549 -.0040102X + (1 . 28 32E-6 )X^ - ( 1 . 3 2 34E-10 ) X^ 

+ (4 . 3908E-15)X‘* 
s 

Q = QO + QKP) + Q2(P)^ 

50 = -400.79 + 1.5801Q - ( 2 . 0 2 54E-4 ) + ( 2 . 4111E-8 ) Q ^ 

-( 8.6737E-13)Q‘^ 

51 = 16.196 -.024333Q + ( 9 . 3484E-6 ) - ( 1 . 2594E-9 ) Q ^ 

+ (4.7522E-14)Q‘* 

52 = -.14758 + (2.359E-4)Q - ( 1 . 0 3 7E- 7 ) + ( 1 . 6 016E-11 ) Q ^ 

-(6 . 3195E-16)Q‘* 

S = SO + S1(R) + S2(R)^ 

MAXIMUM RANGE CRUISE TIME AND SPEED 
AT CONSTANT ALTITUDE PROGRAM 

BO = -1 + (5.0794E-3)H - ( 1 . 39 6 8E- 3 )H^ + (8.254E-5)H^ 

-(1. 2698E-6)H‘^ 
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B1 = .05 + .0015159H + (1.123E-4)H^ - ( 3 . 4921E-6 )H^ 

+ ( 7.9365E-8)H‘^ 

N = BO + BKG) 

BO = .47803 + .0013417D + ( 6 . 2 28 7E-6 ) - ( 1 . 6 2 61E- 8 ) D ^ 

+ (1.6438E-11)D‘' 

B1 = .08217 + (4.1209E-4)D - ( 4 . 557 7E- 6 ) + ( 1 . 6 77 7E-8 ) 

-(2 .OOIE-IDD*^ 

B2 = (4.2143E-4) - ( 9 . 4 39 7E-5 )D + ( 1 . 2 646E-6 ) - ( 4 . 8 5 37E-9 ) 

+ ( 5.7222E-12)D‘* 

B3 = -(6.6767E-4) + (8.4671E-6)D - ( 1 . 0501E- 7 ) + ( 3 . 6 3 8 2E-10 ) 
-( 3.7828E-13)D‘* 

M = BO + BKN) + B2(N)^ + B3(N)^ 

Ml = M -[(60-T) (2) (M)/1200] 

V = (710) (Ml-. 14) + 100 -E 

T1 = Dl/V 

FUEL REQUIRED FOR MAXIMUM RANGE CRUISE 
AT CONSTANT ALTITUDE PROGRAM 

BO = 4.54 -.16444A + .0033932A^ - ( 1 . 02 8 3E-4 ) A^ + ( 1 . 9 2 6E-6 ) A“* 
-(1. 3757E-8)A^ 

B1 = (3.22E-9) -(3.6664E-3)A + ( 8 . 9 3 3 8E-4 ) A^ - ( 5 . 59 39E-5 ) A^ 

+ d.4593E-6)A‘' - ( 1 . 3281E-8 ) A^ 

B2 = (6E-4) + (1.1203E-4)A - ( 2 . 3 3 5 8E-5 ) A^ + ( 1 . 45 36E-6 )A^ 

-( 3. 7144E-8)A‘" + ( 3 . 3334E-10)A^ 

N = BO + BKG) + B2(G)^ 
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BO = -(2. 5399E-3)D + ( 9 . 7 29 9E- 5 ) - ( 2 . 3 516E- 7 ) D ^ 

+ (1.4251E-10)D‘* 

B1 = 2 + (4.2388E-3)D + ( 1 . 2326E-5 )D^ - ( 1 . 0298E-7 ) 

+ (1. 7277E-10)D‘* 

L = BO + Bl(N) 

F = L/V 
R = (F)(T)/60 

MAXIMUM RANGE CLIMB AIRSPEED SCHEDULE 

S = 405.56 -.79075D + .00113820^ - ( 4 . 1018E-7 ) 

M = .86 -(2.1634E-3)D + ( 7 . 6 5 8 2E-5 ) - ( 1 . 1344E-6 ) 

+ (7.2125E-9)D‘* -( 2 . 3035E-11)D^ + ( 3 . 6 5 88E-14 ) D® 

-(2. 3062E-17)D^ 

TAKEOFF AIRSPEED PROGRAM 

U1 = 54.023 + (3.4787E-3)G - ( 1 . 94 75E- 8 ) 

U = U1 + [(26-P)/2] 

VO = -1917.1 + 61.604U -.70348U^ + .0035661U® - ( 6 . 6 5 78E-6 ) U“ 

VI = 76.824 -2.4517U + . 028779U^ - ( 1 . 4 7 5 3E- 4 ) U® + ( 2 . 7872E-7)U‘* 
V2 = -.72239 + .023415U - ( 2 . 79 8E-4 ) U^ + ( 1 . 4 596E- 6 ) U® 

-(2. 807E-9)U“ 

V3 = VO + VKR) + V2(R)^ 

MAXIMUM REFUSAL SPEED PROGRAM 

BO = -43.01 + 6.761G -.35159G^ + .0080545G® - ( 6 . 776 9E- 5 ) G** 

B1 = 26.312 -3.8382G + .20326G^ -.047022G® + (3.994E-5)G^ 
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B2 = -4.9639 + . 72723G -.0387210^ + ( 8.985E-4)G^ - ( 7 . 6 38E- 6 ) G"* 

B3 = .30288 -.044855G -.0023921G^ - ( 5 . 5549E-5 )G^ 

+ (4. 7217E-7)G‘* 

R = BO + Bl(E) +B2(E)^ + B3(E)^ 

BO = -11.412 + 62.185L -9.0037L^ + .64921L^ -.017455L‘^ 

B1 = -.2811 -4.2012L + .70377L^ -. 058693L^ + .0017461L'* 

M = BO + Bl(R) 

OPTIMUM ENDURANCE ALTITUDE PROGRAM 

BO = 55.333 + .073076D - ( 9 . 78 36E-4 )D^ + ( 3 . 5015E-6 )D^ 

-( 3.9782E-9)D“' 

B1 = -1.1 -(8.0597E-3)D + ( 8 . 009 7E-5 )D^ - ( 2 . 8 8 36E- 7 )D^ 

+ ( 3 . 3032E-10)D“ 

B2 = (6.6667E-3) + (1.2541E-4)D - ( 1 . 40 39E-6 ) 

H = BO + BKG) + B2(G)^ 

CRUISE CEILING PROGRAM 

BO = 85.118 -.29117D + .0030434D^ -( 1 . 2851E-5 )D^ + ( 1 . 66 21E-8 )D‘* 

B1 = -2.7877 + .025635D - ( 3 . 306 3E-4 ) + ( 1 . 416 2E-6 ) 

-d.8 34 3E-9)D‘* 

B2 = .063327 - ( 8 . 5 2 89E-4 ) D + ( 1 . 0 814E- 5 )D^ - ( 4 . 6 514E-8 ) 

+ (6 .0606E-11)D‘' 

B3 = -(6.0468E-4) + (9.0826E-6)D - ( 1 . 143E-7 )D^ + (4 . 9304E-10 )D^ 
-(6 .4567E-13)D'" 

H = BO + (Bl)G + (B2)G^ + (B3)G^ 



APPENDIX D 



HP-9830 Programs and Lists of Variables 



I REH THIS PROUPHM CfiLCULfiTES THF FUEL FLOIJ AND LBFUEL 'MhUTIC ftL MILE FOR AH 
R REN A-7E Fl.’f'IHC A LON LEVEL HISSIOH AND IS DEPENDEHT UH 4 VARIABLES — 

3 REN GROSS NEIGHT.DRAG COUMT-NACII HIINBER- AND TENPERATUPE' CENTIGRADE) 

10 PRINT "ENTER GPObS NT . DRAG CT- NACH #. AND TENP-: CENT ) " 

II PRINT 
12 PRINT 

20 INPUT G>D<N,T 
4? G=G/1O00 

5 0 M 1=0, 3 o y 1 5 + 0 , 0 0 4 2 3 0 1 ^ G 
54 GOSlJB SOO 

5d 1=0 

53 G03U8 600 

6 0 S 2 = S 

70 IF Si; 32 THEN lOO 
••^0 S=S2 
■^5 GOTO 300 
100 1=1 
110 G03UB 600 
120 33=3 

130 IF 3K33 THEN 200 

140 32=33 

150 1=1+1 

160 GO SUB 600 

170 GOTO 120 

200 1 1 = ( 3 1 -32 > . ( 33-32 > 

210 M1=M 

220 1=1-1+11 

221 I=INT(I> 

222 G03UB 60O 

223 32=3 

224 1=1+1 

225 G03UB 60 O 



226 33=3 

240 GOTO ?O0 

235 PRINT 

236 PRINT 

300 R=S + 2 + ( 4 . 3732E-O3 + 0. 027743*D> 42 

301 R3=R 

302 R1=2+INT':;R/2> 

304 R2=Rl+2 

306 J=1 

303 IF J=2 THEN 311 

3 0 '? R = R 1 

310 GOTO 319 

311 R=R2 

J 1 9 B U = D . 62 5 — 1 . 9 III F' + j> , y 252'^ F’ v 2 ■" 1 . U 7 6 1 1- R 1 2 + y , 17675 + R t' 4 
320 B0=B0-0. 0 1 3095-^P T 5,f 3 . 526E-04^P T 6 

330 Bl=205. 3012-243. 9317^P+91 . 66355-Pt2- 1 5 . 552 1 3^R +3+ 1 . 224^32-Rt4 
34 ri Bl=Bl-0. 0395333-Pt-5t-2. 896335E-04^P t-6 

350 B2=-1052. 123+1231 . 24t-R-437. 4233 + P 12 + 91 . 6522 + R 1 3-3 . 662962"-R T 4 + 0 . 3953974 + RT 5 
I: b 0 B 2 = E: 2 - O . O O 6 9 O 5 5 5 R T 6 

370 B3=1630. 142-1950. 1 ::9-R + 733. 35 1 3^R T 2- 1 52 . 5733^P T 3+ 1 5 . 03319 + Rt -4 
•i* C' ij E* i' = E' -■ — y ■ t' 2 1 4 1 1* 9 H 1 5 + y . y 1 7 y 7 + R F* 

390 B4 = -364 . 6375+1 OOO . 443 + P-403 , 745 1 f-2 + 30 . 033 1 4:**R T 3-3 . 03953+P T -4 

4 y y B 4 = El 4 + y , 9 31 2 5 2 7 1 R t' 5 — 7 . 7 2 O 6 1 7 E — 0 3 + F! 1 6 
430 ri=B0 + Bl+M + B2 + Nt2+B3^Mt-3 + B4trMt4 

440 IF J=2 THEN 430 
450 N1=N 
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OK IHTEPHEDIFITE « BY LI HEAR IMTERPOLATIOH 



455 J=2 
4E0 GOTO 311 
475 R=2: 

4S0 H2=M 

4'?0 H=H 1 + ■: N2-H 1 ■ - C P3-R 1 ■ 2 

500 REH COHPLETED CflLCULflTlOH 

510 P=4. 9746*H + 7. '?043E-0t>H r-2 



52U 


\\4 = (6 


. 4375+0 


. 010428 + T-6. o925E-0t>TT-2 + 4. 9127E 


5 50 


F = -;'0. 




-1000 


539 


F=INT 


(F> 




540 


PRINT 


"GROSS 


UT=" G^1000 


541 


PRINT 


"TS="S 


"IiC = "D"M="M 


542 


PRINT 


"TEMP= 


"T 


543 


PRINT 


"REF # 


= "R3 


544 


PRINT 


"N="N 




545 


PRINT 


"LBFUEL NM:="P 




PRINT 


"FUEL 


FLOW= "F 


551 


PRINT 






555 


GOTO 


10 




»500 


B0 = 22 


. SI 9-31 


. 734 -1+41 . 33+1 1-2-5. 0953+ n-3 


*410 


Bl=-1 


54.98+2 


17. 5UI-281. 7 5:-I r-2 + 35. 905tl t-3 



620 B2=405 . 00-525 .56+1 +607 . 49+ I r2-yS . 737+1 t3 
630 B3=-445. 62+542.93*1-611.55+142+92. 394+1 t3 
640 B4= 184. 73-204. 42+1 +225. 39+1 1-2- 35. 139+1 f 3 
650 S=B0+B 1 +M 1 +B2+t'1 1 42 +B3 + t'1 1 4.I;+E;4+M 1 44 
660 RETURN 

300 A0=-2. 3237-0 . 263 1 6 + D+0 . 0073327 + D 42-7 . 5 1 3E-05 + D4'3+3 . 5396E-07 +D44- 
310 A0=h0-7. 73E-10+D45+6. 4624E-13-D46 

320 Hi =4. 835+1 . 0956-D-0. 030653+D42+3. 1 9 1 2E-04 + D43- 1 . 5276E-0+I*Ii 44 
330 Al=fll+3. 403E-y9*D4j— 2 . 3692E-12-D46 

340 02 = 10.234-1. 07 19 + D+O. 03 1 094+D4 2-3. 2373E-04 + D4 3+ 1 . 595E-06+D4 4 

350 02=02-3. 6O09E-O9-D45+3 . 0634E-12+D46 

360 3 1 =00*0 1 *H 1 +02+N 1 r-2 

370 RETURN 

330 END 
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List of Variables for Program 1 



Variable 


Definition 


G 


Gross weight (lbs.) 


D 


Drag count 


T 


Temperature (°C) 


M 


Mach number 


Ml 


Result of lower graph, Figure B1 


I 


Guidelines , numbered top to bottom consecutively 


S 


Transfer Scale calculated as function of I 


SI 


Transfer Scale calculated as function of D 


S2 


Transfer Scale calculated for upper guideline 


S3 


Transfer Scale calculated for lower guideline 


11 


Relative Transfer Scale location between guidelines 


R,R3 


Reference number 


R1 


Even reference number below actual reference number 


R2 


Even reference number above actual reference number 


J 


Integer counter 


N 


Result of lower graph, Figure B3 


N1 


Result of lower graph. Figure B3 for R1 


N2 


Result of lower graph. Figure B3 for R2 


N4 


Result of lower graph. Figure B4 


AO ,B0 , 
A1,B1 


Coefficients 


A2,B2, 

B3,B4 


Coefficients 


P 


Pounds of fuel per nautical mile 


F 


Fuel flow 
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Program 2 



1 REM THIS PPluIRmM ChLCULhTES THE Thf EuFF DISTRMCE PEijUIPED FOR hH H-7E 

2 REM IT IS riEPEHUEHT S VhRIhBLES — 

3 PEM GROSS HEIGHT. PHUV HLT ITUDE » TEMP . DRhG COUHT . RELhT I VE HUM ID I T*,M-J I NDS 

4 REM RHUV SLOPE « CE'.NTER OF GRHVITV LOChT I ON » FLOPS « HMD DOUBLE DhTUM SThTUS 

9 PRINT •* INPUT HLT.TEMP* DC.Gir* 

10 INPUT H,B»D.G 

11 1=50 

12 L=10 

13 N=1 

14 p=27 
20 R=25 

100 B0=13. OSb-0. 000171 13^H-2. 0b55E-07-^R T2+3 . bSblE-1 1 ^RtS 

101 B0=BO-2. 41 5bE-15'^R t 4 

1 10 Bl=-0. 045b:; j-7. 393l£-0b'^R+3. 7545E-0’?^RT2 

111 B 1 - B 1 - 9 . 7 0 S ;E: E - 1 2: - R 1 3 + b . 9 9 7 E - 1 7 ^ R 1 4 

120 B2 = -0. 001317-;::. L55SE-07^R^•4 . 0739E- 10>R t'2 

121 B2=B2-S. 54:E:E-14-R r3f5. 49b4E-l'3^RT-4 

1 3 0 £; !• = - 1 . 9 U 9 7 E - U H + 1 . >; b 7 1 E — U ::: •i' R - 9 . 4 b 9 4 E - 1 2 R L + L . U 4 2! 4 E - 1 5 R 1' 3: 

140 B3=B3-l . 4E17E-19^Rt4 

150 C=B0 + Bl-B + B2-^B'T2fB3^Bt3 

IbO IF D=1 THEN 190 

170 E=0. 5417S+0. b5S7b^C 

l:B0 GOTO 200 

190 E= 1 , 9773 + 0 . Jb59:5+C 

200 G0=-4. :E::::9bE t-05+:E:. 4974E + 0 1 ^G-5 . 7:E:5bE-03 + G 1 2+ 1 . 9373E-07^G 1 3- . 1 744E- IZ^G T 4 
210 G0=G0 + 2 . 044 bE- 1 7-G t5 

220 G 1 =5 . Ji'bE 1 E+U4- 1 , 0 1 46E + 0 1 ;^G + 6 . ;:;:E:07E-y 4 t '2- *' , 292E - u:;: +U I :;+ Ei . 73:E:7E- 1 3 +G l' 4 
230 Gl=Gl-2. 39b4E-l:3+G t'5 
240 H=G0+G1+E 
250 J=0 

2bO IF K40 THEN 2t:0 
2 7 0 J = ( 1-40 •• 1 0 0 0 
2:50 K=H+J+H 
2:55 IF L=0 THEN 340 

290 L0=b. 7124E+01 r:?. 95O9E-0 1 ^K + 2 . 330bE-05^Kt2 - 1 . b254E-09-^K T'3 + 3 . 372:E:E- 1 4+F T-4 
300 Ll=-9. 0995-1 . 0:f:5bE-02-K + 2. I 754E-07-^KT2-2 . 5327E- 1 1 -K T-3+ 1 . 1 97E- 1 5 + KT4 
blu L2=l . 47:E:L'E-0 i 1 -2 . 1 bbbE-0tb ••■£ + 5 . 4274E-09 + K + 2-2 . 7:5 1 7E- 1 3 + K I' '3 + 9. 3077E- 1 :5-:^K. 1‘4 
320 M=L0+Ll+L+L2+Lt2 
330 GOTO 350 
340 M=K 

350 :<0 = 4. 5704E + 01+9. 3429E-01-M + 2. 22b5E-05 + Mt2-2. 33:3E-09 + Mt3-r7. 941E-14 + M44 
J LI X 1 = 7 . 9472+1 . 4 '+ 1 4 E — U 2 M + 9 . U 7 U E - Lt b fTt* 2 - 7 . 1 2 2 5 E — 1 U ^ M 3 < + . U b 3 4 E - 1 4 + M t' 4 
370 X2=5. 3blb-S. 513bE-03-M^-3. 59l4E-0b + MT-2-4. 5932E-10 + M r3+l . 9S:59E- 1 4^M t4 
3:50 X = X0+X1+H+X2+HT2 

390 00 = 2. bO42£ + 03-2 . 1 B94^X+ 1 . 09 1 5E-03^L: t 2- 1.111 9E-07 + X t 3 + 3 . bb2E- 1 ^ 4 

400 Ql=-1 . 7573E + 02+2. 260 i E-0 1 -^7:-7 . 5225E-05 + Xt 2 + 7 . 70 1 :E:E-09-^Xt ;5-2 . 5437E- 1 3+X t4 
410 02=2. 0549-4. O102E-O3^7;+l . 2032E-Ob^X t-2- 1 . 3234E- 10^:=:t3 + 4. 390SE- 1 5+X + 4 
4 2 0 0 = Q 0 + 0 1 -T F‘ + ij 2 + P 2 

42L1 2U = — 4 . L<U79E + U2+1 . 5:301+0 — 2. Li254E— L14 + 0’I 2 + 2 . 411 1 E — O::*' — . b737E— 1 3“*^L!'1'4 
440 S 1 = 1 . b 1 9bE + Li 1-2.45 33;E-02^O + 9 . 3404E-0b^Q42- 1 . 2594E-09 + Q t- 5 h 4 . 7522E- 1 4:^-Q t-4 
450 S2 = -l . 4753E-01 f-2 . 359E-04 +0- 1 . 037E-07-U 42+ 1 . BulbE-l 1 3-6 . 3 1 95E- 1 b + Qt4 

460 S = S0 + SUR+S2 + RT'2 
470 S=INT(S> 

479 PRINT “ FOR" 

400 PRINT ••GN = -’G" RLT=“R” TEMP = "B" DC = " D ” PH= ” T' HDND= “L 

402 PRINT *’RNNV SLP = “N'''. CEN GRRV= " P ” FLRPS= “ R 

403 PRINT 

530 PRINT "TRKEOFF ROLL DI0T='*S 

531 GOTO 9 

532 END 



74 



List of Variables for Program 2 



Variable 


Definition 


A 


Runway Altitude (feet) 


B 


Temperature (°C) 


D 


Double datum status (1 indicates "with") 


G 


Gross weight (lbs.) 


I 


Relative humidity (%) 


L 


Headwind (kts.) 


N 


Runway slope (%) 


P 


Center of gravity (%) 


R 


Flap position (degrees) 


C 


Result of upper graph, Figure B5 


E 


Takeoff factor 


H 


Unadjusted ground roll distance, Figure B6 


J 


Adjustment factor due to relative humidity 


K 


Ground roll distance (GRD) adjusted for relative 
humidity 


M 


GRD adjusted for wind 


X 


GRD adjusted for runway slope 


Q 


GRD adjusted for the center of gravity location 


s 


True GRD (also adjusted for flap position) 


BO , GO ) LO ) 
X0,Q0 


Coefficients 


B1,G1,L1, 
XI, Q1 


Coefficients 


B2,G2,L2, 

X2,Q2 


Coefficients 


B3,S0, 

S1,S2 


Coefficients 
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Program 3 



1 REM THI'i PROGRhM ChLCULhTES THE ft-7E Mfl.'lIHUM PfttIGE' fi I P.SPEED HlID 

2 REM TIME UF FLIGHT AMD li DEPEMDEHT UM 6 VflRlHBLE':: -- 

:3 REM GROSS WE I GHT - ALT I TUDE « DRAG COUHI - TEMPERATURE W I HDS AND D I STANCE 
S PR I NT "I HPUT GW . ALT , DC « TEMP ^^C > ! HDWD • D I STANCE " 

10 INPUT G. H^D.'TjL.DI 
3y G=G.-" 1 000 
40 H=H/1000 

50 A0^-l+5. 07'?4E“0:::^H“1 . 3‘?G8E-03±H . 49 1 2E-0G"rHt3 ^-7 . 33G5E-03^Ht4 

E0 R1 =0. 05 + 0 . 00 15 159-^ Hi- 1 . 1 23E-04-rHt'2- 3' . 4921 E-0Ei^Ht 3 + 7 . 93o5E~03irHt4 
70 H = H0 + R1-:G 

30 B0=0 . 47303-0 . 00 1 34 1 7-^D + E . 2237E-0E-^D 1 2- 1 . 323 1 E-03 + D-t 3-r i . 343$E- 1 1 -^D t 4 
35 B1=0. 03217+4. 1209E'04iTD-4. 5577E-0t.-Dt2+l . 3777E-03 ^Dt 3-2 . 03 1 E- 1 1 -D t4 

9 U B 2 = 4 . 2 1 4 ji E - 0 4 - 9 . 4 cJ 9 7 E — 0 5 D + 1 . 2 3 4 3 E — 0 3 D 4' 2 - 4 « IJ 5 -3 7 E - 0 9 D 4’ 2 + 5 . 7 2 2 2 E - 1 2 + D 1‘ 4 
95 B3=-3 . 3737E-04 + 3 . 437 1 E-03^D- 1 . 050 1 E-07 + D4-2 + 3 . 3 332E- 10 + D t-3-3 . 7323E- 13 + D4 4 

1 0 0 M = E! 0 + B 1 + f J + B 2 ^ N 1 2 T B 3 ir N r- 

1 10 ( '.:30-T)+2 + M. 10+120> 

1 20 V = 7 1 0+ M-0 . 14 :■ + 1 00-L 
130 T1=D1.-V 
135 V=IHT‘.V> 

140 PRINT 

150 PRINT •■FOR" 

130 PRINT "GN=”G'‘ RLT = “H“ Di: = ’• D '• TEMP = ‘‘ T“ HDHD= “ L " D 1 3T= " H 1 

131 PRINT 

132 PRINT GROUND •5REED=’‘ V’* TIt4E OF FLIGHT = “T1 

170 END 
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List of Variables for Program 3 



Variable 

G 

H 

D 

T 

L 

D1 

N 

M 

V 

T1 
AO, BO 
A1,B1 
B3 



Definition 

Gross weight (lbs.) 

Altitude (ft.) 

Drag count 
Temperature (°C) 

Headwind (kts.) 

Distance to fly 

Result of first chart, Figure B8 

Cruise Mach number (adjusted and unadjusted for T) 
Ground speed (kts.) 

Time of flight 
Coefficients 
Coefficients 
Coefficient 
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Program 4 



10 REM THIS PRijCRhM CRLCULhTES FUEL REuUIRED FOP MR’' PRNGE RT CONSTRHT 

11 PEM RLTITUDE FOR RM R-7E RND IS riEPEHDENT ON S VRPIRBLES -- 

12 PEM GROSS HEIGHT. RL T I TUDEiDPRG COUNT ^ TRUE RIPSPEED. RND T I ME < M I NOTES > 

20 PRINT “ENTER GROSS UT.RLT^DRRG CT ^ TR3 ^ T I ME < M I NOTES > “ 

21 PRINT 

22 PRINT 

S0 INPUT G.R.D^ViT 
S.5 PRINT “GROSS HT = "G 
3b PRINT ■'RLTITUriE=“R 

37 PRINT “IiRRG i::OUNT = ”D 

38 PRINT “TRUE RIPSPEED=“V 

39 PRINT “TIME OF FLIGHT=“T 

40 G=G/1000 
50 R=R-'’'1000 

bO B0=4. 54-0 . 1S444 tR + 0, 0033932 -^R T ^- 1 . u283E-04-^Rt3+ 1 . 92bE-06^R T 4- 1 . 3757E-08 -^RtS 
70 Bl=-3. bbb4E~03^R + 8. 9338E-04 + R T 2-5 . 5939E-05^R T '3+ 1 . 4593E-0b^R T-4- 1 . 328 lE-08;^R t'5 
80 B2=bE-04 + l . 1203E-04.-R-2. 3358E-05^RT2+1 . 453bE-06'^R 1 3-3 . 7 1 44E-08^Rt4 
85 B2=B2 + 3. 3334E-10-RT-5 
90 N=B0+B1^G+B2-Gt'2 

100 R0 = -2. 5599E-03^D+9. 7299E-05-Ii 1 2-2 . 35 1 bE-07^Ii T 3+ 1 . 425 1 E- 1 0*ri T-4 
110 Rl=2+4 . 2388E-03^Ii+ 1 . 232bE-05"-Ii T' 2- 1 . 0298E-07-D t 3+ 1 . 7277E- 10^0 ^4 
120 L=R0+R1+N 
130 F=L^V 

140 R = F ^ T b 0 

141 L= '• I NT < 1 000 ■' / 1 000 

142 F=INT(F> 

143 R=INT(R> 

147 PRINT 

148 PRINT 

150 PRINT “LBFUEL. NM=“L“FUEL FL0N=“F 

155 PRINT “FUEL REOU I R.EIi=“ R 

15b PRINT 

157 PRINT 

IbO GOTO 20 

170 END 
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List of Variables for Program 4 



Variable 

G 

A 

D 

V 

T 

N 

L 

F 

R 

BO ,A0 
B1,A1 
B2 



Definition 

Gross weight (lbs.) 

Altitude (ft.) 

Drag count 

True airspeed (kts.) 

Time of flight (minutes) 

Result of first chart, Figure B9 

Pounds of fuel per nautical mile 

Fuel flow 

Fuel required 

Coefficients 

Coefficients 

Coefficient 
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Program 5 cr ^ 

J 



1 PEM THIS PROGRhM ChLCULhTES THE CLIMB HlRSPEED OF hH H-?£ 

2 REM • INDICHTEn RIR’SPEED BELOH EO^OOOV' 

3 R E M M H C H N U M B E R H B 0 V E 2 0 ? 0 0 0 ’ ) 

10 n=0 

12 PRINT "CLIMB HlRSPEED SCHEDULE" 

15 PRINT "DRhG CT climb hIRSPEED CLIMB MhCH" 

1 6 P R I N T " I H S T ij 2 O U O O ^ > ( R B ij E 2 O O O U • " 

20 S=403 . SE-0 . 790P j->-D + 0 . OO 1 1 3S2*D t 2-4 . 101 BE-07^DT'3 

21 S=INT-:S..' 

30 M=0. 3E-2. 1B 34E-03 <^0 + 7. SS82E-’0S^DT2- 1 . 1 344E-0r.:-D t 3 + 7 , 2 1 2 jE-09-D 1 4-2 . 3 

4 U M = f1 + 3 • B j 3 S E ” 1 4 ^ D l' B - 2 . 3 0 6 2 E - 1 7 ^ D t* 

42 M=M.fl0O0 

44 M=INT<M> 

4B M=M-- 1000 
55 PRINT DiSmM 
B0 D=D+30 

70 IF D<310 THEN 20 
80 END 



35E-1UD T 
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List of Variables for Program 5 



Variable 


Definition 


D 


Drag count 


M 


Mach number 


S 


Calibrated airspeed (kts.) 
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1C ? M 
•tE 



13 



400 

401 

402 

493 

499 

500 

501 

502 

503 

504 

530 

540 

550 

560 

570 

5 y y 

590 

600 

610 

620 

625 

630 

63 1 

635 

640 

650 

651 

652 

655 

656 

657 

660 

662 

663 

669 

670 

710 



Program 6 



REM THIS PROGRhM ChLCULhTES THE TRKEOFF RIRSPEED OF RN R-7E 

REM UNDER VRRVIHG GROSS HEIGHTS^ FLRP POSIT IONS 7 

REM RND CENTER OF GRRVITV LOCRTIONS 

R = 20 

P = 20 

G=20000 

PRINT ’’FOR GROSS NEIGHT="G 

PRINT 

PRINT 

PRINT "FLRPS CG TRKEOFF RIRSPEED’* 

Ul=5. 4023E+01+3. 4737E-03+G-1 . 9475E-03-G t-2 
U=Ul + ':26-P>.-2 

V M = — 1 , 9 1 7 1 E "i" U 6 . 1 6 0 4 E U 1 U *” 7 , y S 4 iri E ~ U 1 t J 2 + . j 6 6 1 E “ U IJ J ~ 6 . 6 5* 7 E *” y 6 Li 4 

VI =7 . 6324E + 01-2. 4517-U + 2. 3779E-02-lJt-2- 1 . 4753E-04:^U T 3+2. 7872E-u7-rU T 4 

V2 = -7. 2239E-01+2. 34 15E-02 + U-2. 793E-04-Ut2-^ 1 . 4596E-06+Ut '3-2 . 307E-O9 + U f 4 

V 3 = 0 + 1 ^ R + V 2 + R T- 2 
V4=INT(V3) 

PRINT RiP?V4 
R = R + 5 

IF R;:-40 then 6 2:0 
GOTO 530 
P=P + 3 
R=20 

IF P>35 THEM 650 

GOTO 530 

G=G+300O 

R=20 

P=20 

PRINT 

PRINT 

PRINT 

PRINT ••FUR GROSS NEIGHT = "G 

PR I NT 

PRINT 

IF G >42000 THEN 710 

GOTO 530 

END 
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List of Variables for Program 6 



Variable 


Definition 


R 


Flap position (degrees) 


P 


Center of gravity (%) 


G 


Gross weight (lbs.) 


U1 


Unadjusted takeoff airspeed 


U 


Takeoff airspeed adjusted for center of gravity 


V4 


Actual takeoff airspeed (adjusted for flap position) 


VO, VI, 
V2,V3 


Coefficients 
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Program 7 



1 REM THIS PROi'.RmM CHLCULhTES THE MhI'IMUH REFUSfiL SPEED 

2 REM FOR hM ft-?E USING ANT I -SKID 

3 REM IT IS DEPENDENT ON 5 VARIABLES — 

4 REM GROSS HEIGHT- TEMP j RNMY LENGTH -RNHV ALTITUDE- AND DOUBLE DATUM 

9 PRINT "INPUT ALT- TEMP- RNHY LTH- GW> DOUBLE DATUM" 

10 INPUT A-B-L-G-D 
1j G=G-1O0O 

20 L=L''lOOO 

70 PRINT "ALTITIJDE="A 

71 PRINT "TEMP="B 

72 PRINT "RNWY LTH="L*1OO0 

73 PRINT "GROSS NT = "G-^1OO0 
75 PRINT "DD="D 



100 


B0 = 1 3 


. 0 


36^0 , 000 


1711 


3^R-2. 06 


55E 


-07 + RT- 


2 + 2i . 6 3 6 


lE-1 


l*At'3 




101 


B0=B0 


“ 2 


. 4156E-1 




4 














110 


Bl=-0 


. 0 


45635-7. 


8931 


E-06+R+3 


. 1 J 


45E-09 


-rRt-2 








1 1 1 


B1=B1 


”9 


.7080E-1 


3^Ht 


T: f 6 . 9 9 7 E 


-17 


T^Rt4 










120 


B2 = -0 


. 001317-3. 


2558 


E-07i^R + 4 


. 07 


39E-10 


+ RT2 








121 


B2=B2 


-0 


.548E-14 


+ Rt3 


+5. 4964E 


-18 


-RT-4 










r:-:0 


B3 = -l 




097E-05+ 


1 . Jb 


71E-08+R 


— 9 . 


4694E- 


12^RT2+ 


2. 04 


34E-1 


5-^RT-3 


140 


B3 = B3 


-1 


.4617E-1 


9 + Rt 


4 














150 


0 = B0 + 


B1 


-B+B2-BT 


2+B3 


+ BT-3 














160 


IF D = 


1 


THEN 190 


















170 


E=0.5 


41 


78+0. 658 


76 +C 
















180 


GOTO 


200 


















190 


E = 1.9 


r 1 


Ti+0 . 5659 


8+C 
















200 


B0 = -4 


4' • 


01+6 , 76 1 


^G-0 


. 35159 ^^G 


1 2 + 


0 . 0 0 8 0 


545+GT3 


-6 . 7 


763E- 


05-Gt4 


2 1 0 


Bl=26 


, '3 


12-2. 838 


2^G + 


0 . 20326+ 


G 12 


-0 . 004 


7022+G t- 


:5 + 3 . 


994E- 


05-^G'T-4 


220 


B2=-4 


, 0 


639+0 . 72 


i kZ! 'S ^ 


G-0 . 0387 


21 - 


GT2+8. 


985E-04 


--G t-3 


— 7 . 63 


8E-06+Gt4 


230 


B3=0. 


30 


288-0. 044855 


+ 1 j + 0 . 0 0 2 


3 ■? 2 


1 :h: t'2- 


5.5549E 


-05 + 


Gt-3 + 4 


. 7217E-07^GT4 


240 


R=B0 + 


B1 


^E+B2^Et 


2 + B 3 


+ET3 














250 


B0 = -1 


1 . 


412 + 62. 1 


85-L 


-9 . 0037-^ 


Lt2 


+0 . 649 


21+Lt3- 


0 . 0 1 


7455 + 


Lt4 


260 


Bl=-0 


• 2 


811-4.20 


12-L 


+0 . 70377 


+ LT' 


2 - 0 . 0 5 


8693^Lt 


3 + 0 . 


00174 


61-Lt4 


270 


M=B0+ 


B1 


+ R 


















271 


M=INT 


ai 


> 



















272 PRINT 

275 PRINT ■■ TAKEOFF FACTOR ="E 
27F. PRINT 

230 PRINT "MAX REFUSAL SPEED = "«M 

231 PRINT 
230 GOTO 3 
300 END 



1 



•-::ThTUS 
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List of Variables for Program 7 



Variable 

A 

B 

L 

G 

D 

C 

E 

R 

M 

B2,B3 



Definition 

Runway Altitude (ft.) 

Temperature (°C) 

Runway length (ft.) 

Gross weight (lbs.) 

Double datum status (1 indicates "with") 
Result of upper chart, Figure B5 
Takeoff factor 

Result of first chart, Figure B12 
Maximum refusal speed (kts.) 

Coefficients 
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Program 8 



1 REM THIS PROGRRM CRLCULRTES THE OPTIMUM ENDURRNCE RLTITUDE 

2 REM OF RN R-7E RT VRPV I MG GROSS HEIGHTS RND DRRG COUNTS 

4 DIM BCS3 

5 G=IS 

6 D=0 

10 PRINT ■‘OPTIMUM ENDURRNCE RLT ** 

20 PRINT •■ GROSS HT DRRG CT OPT END RLT" 

50 G=G+3 

SO BC 3 ]=55 . 333 + 0 . 07307S ^D-9 . 7S3SE-04^D t-2 + 3 . 50 1 5E-0G+D +3-3 . 37S2E-09 + D t-4 

90 BC I ]=-l . I -3. 0597E-03^D + 3 . 0097E -05-^D T-2-2 . 383SE -07 t 3 + 3 . 3032E-10+D + 4 

100 BC 2 ]=S. 6667E-03+1 . 254 1 E-04^D“ 1 . 4039E“0S^£it2 + 5 . 2032E-09^Dt3-G . 0213E-12^D + 4 

110 H = BC3] + BC 1 ]^G + BC 2 ]^Gr-2 

115 Z=INT<'H+iOOO) 

IIS N=G+1000 

119 PRINT XjDkZ 

120 D=D+30 

121 IF B<310 THEN SO 

122 D=0 

123 IF G<45 THEN 50 
140 END 
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List of Variables for Program 8 



Variable 


Definition 


G 


Gross weight (lbs. times 1000) 


D 


Drag count 


H 


Optimum endurance altitude (ft.) 


Z 


Optimum endurance altitude ( integer format) 


X 


Gross weight (lbs.) 


31,B2 ,B3 


Coefficients 
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Program 9 



1 REM THIS PROGRmM CRLCULhTES THE CRUISE CEILING OF- RN R-7E 
S REM UNDER VRRVING GROSS HEIGHTS RND DRRG COUNTS 

4 DIM BC4] 

5 G=19 
R. D=0 



10 


PRINT ■ 


•CRUISE 


CEII 


20 


PRINT * 


•GROSS NT 


50 


G=G+3 






SO 


B[4]=S5. 118-0. 


2-^1 


90 


BC 1 ]=-2. 7877 hO 


. 02‘ 


1 00 


B[ 2 ]=0. UE8327 


— C' > 


105 


B[3] = - 


" b . 04b 3 E 


-04 


110 


H=BC4 


] + &[ 1 ]^G 




115 


Z=INT' 


• H-^IOOO) 




118 


x=G^ii; 


j 0 0 




119 


PRINT 


71 .1 D 1 Z 




120 


D=D+:]:i? 




121 


IF dc: 


51 Cl THEN 


80 


122 


D=0 






123 


IF G<‘ 


45 THEH 


50 


140 


END 







DRRL CT 



CRUISE CEILING** 
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List of Variables for Program 9 



Variable 


Definition 


G 


Gross weight (lbs. times 1000) 


D 


Drag count 


H 


Cruise ceiling (ft.) 


Z 


Cruise ceiling (integer format) 


X 


Gross weight (lbs.) 


B3,B4 


Coefficients 



89 



APPENDIX E 



TI-59 Programs and User Information 



USER INFORMATION FOR PROGRAM 1 
Program: Low Level Cruise Performance 

Number of Steps: 1386 

Computation Time: 90-110 seconds 



STEP 

1 

2 

3 

4 

5 

6 

7 

8 
9 
0 



ENTER 

gross weight (lbs.) 
drag count 
mach number 
temperature (°C) 

read in cards 3 S 4 
drag count 
mach number 
temperature (°C) 



PRESS KEY 
A 
C 
D 
E 

R/S 

C 

D 

E 

R/S 



DISPLAY 

gross weight/1000 
drag count 
mach number 
Transfer Scale 
Unusable number 

Transfer Scale 
mach number 
lb . fuel/nautical mile 
fuel flow 
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000 


76 LBL 


050 


1 j 3 


y 


1 1 j 1 j 


43 


RCL 


0 0 1 


1 1 R 


05 1 


01 




101 


02 


02 


002 


55 


052 


65 




102 


45 


Y ' 


003 


0 1 1 


053 


43 


RCL 


1 0 :j 


03 


J 


004 


00 0 


054 


00 


00 


104 


-7tr 
i J 


- 


005 


Ijlj IJ 


055 


95 


= 


105 


1 j 3 


.ji 


nri6 


ri ij ij 


056 


42 


3 TO 


106 


9 3 




007 


95 = 


057 


05 


05 


107 


05 


CJ 


003 


42 3TG 


053 


“i cr 

w J 


CLR 


103 


0 3 


:3 


rifi9 


00 00 


059 


75 


- 


109 


1 j 9 


'P 


010 


91 R/3 


060 


53 


f’ 


110 


06 


6 


01 1 


76 LBL 


06 1 


02 


4 


1 1 1 


• J 4 


£E 


1 j 1 2 


12 B 


062 


9 "3 


. 


112 


94 


+ / - 


0 1 3 


55 - 


063 


1 j 3i 


3 


113 


07 


T* 


0 1 4 


01 1 


0 6 4 


02 


4 


114 


65 


•4 


0 1 5 


00 0 


065 


ij3 


•J 


115 


43 


RCL 


0 i 6 


1 j 1 j 1 j 


1 j 6 6 


07 


r 


116 


02 


02 


017 


0 Ij 0 


ij b 7 


35 




117 


45 


' , ' . 


0 1 


95 = 


063 


j 3 




1 1 3 


04 


4 


ij 1 9 


42 bin 


069 


02 


4 


119 


35 


+ 


020 


0 1 0 1 


070 


06 


6 


120 


07 


1 


021 


91 R.--3 


071 


03 


.J 


121 


9 3 


, 


022 


76 LBL 


072 


01 


1 


1 22 


07 




023 


1 3 1 j 


073 


1 j 6 


6 


1 4 :! 


1 j 3i 


y 


024 


42 3T0 


074 


65 




124 


52 


EE 


025 


02 02 


075 


43 


RCL 


1 ■ s 


94 


+ - 


02 b 


91 R/3 


076 


02 


02 


126 


01 


4 

i 


027 


76 LBL 


077 


-7 C- 

r 


- 


147 


'j 1 j 


0 


023 


14 D 


1 j 7 Ij 


93 


, 


4 ‘•i i“i 

i 


65 




029 


42 3T0 


Li 7 9‘ 


00 


Ij 


129 


43 


RCL 


flXfi 


ij>! 0;’-: 


ij Ij 1 j 


Oij 


Ij 


130 


02 


02 


03 1 


R/S 


03 1 


07 


i 


1 3i 1 


45 


yx 


032 


76 LBL 


0 3 2 


03 


•J 


132 


05 


”1 


033 




r i i j j 


0 3i 


•T] 


1 3 3i 


■7 *=’ 




i -J L. 






'• 




’ ■* 




034 


42 3T0 


034 


02 


4 


134 


06 


0 


035 


I J 4 iJ 4 


035 


07 




135 


9 3i 


, 


036 


5 •?' ( 


036 


65 




136 


04 


d 


037 


9 3 


037 


43 


RCL 


1 'T! "7 


06 


hi 


033 


1 j 3 3 


033 


02 


02 


133 


02 


4 


039 


1 j 3 Ij 


039 


•J 3 




i 3 9 


04 


4 


040 


1 j y :5 


i j ■? 0 




-i- 


140 


tr 0 


E£ 


04 1 


01 i 


091 


07 


■J7 


14 1 


94 


- 


042 


0 3 3 


092 


93 


, 


142 


0 1 


1 


043 


35 T* 


0 9 ■ j 


05 


5 


14:3 


1 j 3 


■J 


044 


9 3 


094 


01 


1 


144 


65 




045 


0 0 0 


095 


1 j 3 


L' 


145 


43 


RCL 


0 4 6 


1 j 0 0 


096 


D j1 


c c 


146 


02 


02 


047 


04 4 


1 J ’ J 7 


94 


-h ■■ - 


147 


45 




043 


02 2 


093 


05 


5 


143 


Mb 


6 


049 


1 j 9 9 


099 


65 




149 


54 





91 



1 50 


35 + 


200 


ij6 6 


250 


04 4 


151 


53 ( 


201 


52 EE 


251 


75 - 


152 


04 4 


202 


94 


•”i CJ -“| 


0 1 1 


153 


9 3 


203 


Ij 6 6 


C" 

tl. • J 


9 25 


154 


1 J 3 3 


204 


65 X 


254 


1 j 1 j 0 


155 


0 3 3 


205 


43 RCL 


255 


07 7 


15b 


05 5 


206 


0 2 1 j 2 


256 


01 1 


15? 


35 + 


207 


45 yx 


•i cr r» 
i_ r 


0 9 9 


153 


01 1 


20 3 


04 4 


-1 C .-1 

^ ■ J c* 


6 5 X 


159 


9 3 


209 


35 + 


259 


43 RCL 


IbU 


ij 0 0 


2 1 0 


ij3 3 


2 6 1 j 


! j 2 0 2 


Ibl 


1 j 9 9 


21 1 


9 2'> • 


261 


35 + 


162 


05 5 


212 


04 4 


2 b 2 


9 3 . 


163 


! j 6 6 


213 


0 1 j 1 j 


2 6 3 


0 1 j 0 


164 


65 X 


214 


03 3 


26 4 


0 3 25 


165 


43 RCL 


215 


52 EE 


265 


01 1 


166 


1 j 2 0 2 


216 


94 
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USER INFORMATION FOR PROGRAM 2 



Program: Takeoff Ground Roll Distance 

Number of Steps: 1385 

Computation Time: 48-50 seconds 



STEP ENTER 

1 gross weight (lbs.) 

2 pressure altitude (ft.) 

3 temperature (°C) 

4 headwind (kts.) 

5 relative humidity (%) 

6 

7 read in cards 3 S 4 

3 

9 flap position (degrees) 

0 center of gravity (%) 

1 runway slope (%) 



PRESS KEY 
A 
B 
E 

2nd,D 
2nd,E 
R/S 
R/S 
2nd, A 
2nd,B 
2nd,C 



DISPLAY 
gross weight 
Pressure altitude 
temperature 
headwind 
unusable number 
unusable number 

unusable number 

flap position 
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takeoff ground roll 
distance 
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USER INFORMATION FOR PROGRAM 3 



Program: Maximum Range Time and Speed at Constant Altitude 

Number of Steps: . 547 
Computation Time: 20 seconds 
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USER INFORMATION FOR PROGRAM 4 



Program: Maximum Range Fuel Required at Constant Altitude 

Number of Steps: 439 

Computation Time: 16 seconds 
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USER INFORMATION FOR PROGRAM 5 



Program; Military Climb Schedule 
Number of Steps: ■ 195 
Computation Time: 3 seconds 

STEP ENTER PRESS KEY 

1 gross weight (lbs.) A 

2 altitude (ft.) B 

3 drag count C 

4 temperature (°C) E 

5 — - R/S 



DISPLAY 

gross weight/1000 

altitude 

drag count 

indicated airspeed 
(kts. below 20000 ft.) 

mach number (above 
20000 ft.) 
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USER INFORMATION FOR PROGRAM 6 



Program: Takeoff Airspeed 

Number of Steps: 265 

Computation Time: 9 seconds 



STEP 


ENTER 


PRESS KEY 


DISPLAY 


L 


gross weight (lbs.) 


A 


gross weight 


2 


flap position (degrees) 


2nd, A 


flap position 


3 


center of gravity (%) 


2nd,B 


takeoff airspeed (kts) 
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USER INFORMATION FOR PROGRAM 7 



Program: Maximum Refusal Speed 

Number of Steps: 716 

Computation Time: 26 seconds 



STEP ENTER 

1 gross weight (lbs.) 

2 runway altitude (ft.) 

3 temperature (°C) 

4 runway length (ft.) 



PRESS KEY 
A 
B 
E 

2nd,C 



DISPLAY 

gross weight/1000 
runway altitude 
temperature 
maximum refusal speed 
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USER INFORMATION FOR PROGRAM 8 



Program: Optimum Endurance Altitude 

Number of Steps: 229 

Computation Time: 8 seconds 

STEP ENTER PRESS KEY 

1 gross weight (lbs.) A 

2 drag count C 



DISPLAY 

gross weight/ 10 00 

optimum endurance 
altitude (ft.) 
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USER INFORMATION FOR PROGRAM 9 



Program: Cruise Ceiling 

Number of Steps; 299 
Computation Time: 11 seconds 

STEP ENTER 

1 gross weight (lbs.) 

2 drag count 



PRESS KEY 
A 
C 



DISPLAY 

gross weight/1000 
cruise ceiling 
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APPENDIX F 



Microprocessor Program 

This computer printout was reproduced on an IBM-360/67 
Computer. This listing was entered by hand on an INTEL Corpora- 
tion MDS-800 System, stored on a diskette, and sent to the 
IBM-360/67 via a high speed line. 
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